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Abstract 

Experimental observations of ignition in premixed gaseous 
reactants indicate that perfectly homogeneous initiation is 
practically unrealizable. Instead, combustion first sets in, as a 
rule, at small, discrete sites where inherent inhomogeneities cause 
chemical activity to proceed preferentially and lead to localized 
explosions. Combustion waves propagating away from these "hot 
spots" or "reaction centers" eventually envelope the remaining 
bul k . 

This study examines the spatial structure and temporal evolution 
of a hot spot for a model involving Arrhenius kinetics. The hot 
spot, characterized by peaks in pressure and temperature with 
little diminution in local density, is shown to have one of two 
possible self-similar structures. The analysis employs a 
combination of asymptotics and numerics, and terminates when 
pressure and temperature in the explosion have peaked. 
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1. Introduction 

This paper descibes mathematically the birth and growth o-f a hot 
spot, or localized thermal explosion, in a premixed reactive gas. 

Experimental observations on the initiation o-f combustion in 
gases at high temperatures, in shock tubes and elsewhere, have 
demonstrated conclusively that spatially homogeneous combustion is 
essentially an unattainable ideal. In -fact, ignition -first sets in 
locally, in small volume elements at discrete sites, where chemical 
reaction proceeds preferential ly due to inherent imperfections in 
the system. In due course, combustion waves originating from 
localized explosions occuring at these "hot spots", "reaction 
centers" or "exothermic centers" envelope the entire reacting mass. 
The role played by these sites as precursors of more dramatic 
combustion phenomena is revealed with unsurpassed clarity in Urtiew 
and Oppenheim's Cl 3 photographic records of def lagration-to- 
detonation transition in a Hydrogen— Oxygen mixture confined to a 
tube. These photographs show that as the deflagration travels down 
the tube, it accelerates and evolves into a highly folded turbulent 
flame, preceded by a so-called precursor shock. Eventually, an 
exothermic center is formed in the vicinity of the flame, near the 
tube wall. The localized explosion in this center creates a blast 
wave which propagates through the preconditioned mixture behind the 
precursor shock and ultimately evolves into a fully-developed 
detonation. The same feature appears in other modes of detonation- 
initiation, as well as in other geometric configurations. 

The early analyses of reaction-center dynamics are due to Zajac 
and Oppenheim C23 and Meyer and Oppenheim C33. In these studies 
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the reaction center is assumed to be a spatially homogeneous source 
of chemical energy, capable of expansion and separated -from its 
inert surroundings by an impermeable barrier, across which only 
momentum transfer can occur. Either by prescribing a specific 
reaction scheme, or by specifying an energy release profile within 
the center, the above authors were able to compute the resulting 
pressure pulse. 

In this paper the reaction center is treated as part and 
parcel of the reacting medium rather than an isolated entity in an 
inert atmosphere, and is found to have a definite spatial 
structure. The aim of this paper is to describe this structure and 
to study its temporal evolution in a plane, one-dimensional 
framework, under the asssumption that the reactive gas undergoes a 
single, one-step, first-order, irreversible chemical reaction of 
the Arrhenius type. One may argue that the simple overall kinetic 
scheme adopted here is too idealized to be realistic. However, for 
large activation energies, the kinetics does capture an essential 
attribute of most combustion systems, namely, a reaction rate which 
accelerates rapidly with increase in temperature. Thus the model 
is quite appropriate for studying problems, such as the one at 
hand, which owe their genesis to the interaction between 
gasdynamics and chemical heat release at highly temperature- 
sensitive rates. 

The configuration of the system is so prescribed as to provoke 
the development of a single hot spot, and this can be accomplished 
in a variety of ways. For example, the shock— induced thermal- 
runaway studies of Clarke and Cant [43 and Jackson and Kapila C53 
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considered a semi— in-finite expanse of gas ignited by a piston- 
driven shock, thereby creating a hot spot at the piston face. 
Instead, the present work assumes that the gas is confined between 
two parallel planes, and that its initial state possesses a slight 
spatial nonuniformity. (In a practical situation these non- 
uniformities may be caused by a variety of factors, such as 
turbulence, interacting pressure waves, or, in the case of 
condensed explosives, material imperfections.) The mathematical 
model leads to an initial -boundary value problem for the equations 
of reactive gasdynamics. An asymptotic solution is developed in 
the limit of large activation energy, and the analysis is carried 
as far as the end of the localized explosion within the center. 

The subsequent expansion of the center, and the eventual generation 
of a blast wave, will be the subject of a future publication. 

The temporal evolution of the explosion occurs in two stages, 
beginning with the induction stage. Here the state of the gas is a 
small perturbation of the initial state and the underlying physical 
processes are those of linearized acoustics coupled to a weak but 
nonlinear chemical reaction. The reduced equations require a 
numerical solution (see C41 and C5D which exhibits local thermal 
runaway. Induction is followed by the explosion stage, which 
consists of several distinct spatial zones. There is the 
practically frozen outer zone, and a rapidly shrinking inner zone 
or layer in which intense chemical activity leads to an explosive 
growth of temperature and pressure. Nonlinear chemistry is again 
coupled to linearized gasdynamics, but now the linearization is 
about an atmosphere undergoing a spatially homogeneous thermal 
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explosion. As the layer shrinks, it recedes a way from the outer 
zone, thereby creating an intermediate zone which is -frozen in 
time. Although highly nonlinear, the explosion stage is amenable 
to analysis because gasdynamics is o-f secondary importance; 
temporal variations are much too rapid -for the gas to undergo 
significant expansion. 

For the specific reaction scheme under consideration it is found 
that the reaction center can have one of two posible spatial 
structures, depending upon whether the temperature profile within 
the hot spot has a sharp peak or a rounded peak (Figure 1). The 
former typifies hot spots originating at boundaries (e.g., a piston 
face), and the latter those occur ing in the interior of the vessel. 
These structures, which will be referred to as the "Typs B" 
(boundary— type) or "Type I" (internal type) , are both self —similar . 
The former is described below in detail, with only the results for 
the latter given in section 6. In addition to these two structures 
there exists a third, described briefly in the Appendix; it is 
singular and corresponds to very special initial conditions. 

The specific configuration under study here was also examined, 
with similar methods, by Poland and Kassoy E63. Their analysis 
differs from ours in one crucial respect; they considered the 
distinguished limit in which the spatially homogeneous induction 
time at the initial state and the conduction time across the vessel 
are of the same order, i.e., the Frank-Kamanetski i number 8 is of 
order unity, albeit supercri tical . In our analysis the induction 
time is comparable to the acoustic time across the vessel, i.e., 6 
is very large. In physical terms, the explosive mixture being 
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considered here has a -faster reaction rate. 

2. The Basic Equations, and Bstup 

The enuat i ons o*f reactive nasrivnumi r« -for- nlme. nne-Hi m»n*i nnal , 

• • — » v ■» 

unsteady motion are C73 


(2. la) 

P t + up x + 

Pu x - 0, 

(2. lb) 

p(u^ + UU X ) 

+ (i/y>p x = 0, 

(2. lc) 

p<T t + uT x ) 

- c <y-n/y3 <p t + up x ) 

(2. Id) 

p(Y t + uY x ) 

= -w. 

(2. le) 

p - PT, 



where 

(2. it) w * Cl/<j80)3pYexp<0-0/T). 

Here p, p, T, u and Y are, respectively, the gas pressure, density, 
temperature, velocity and reactant mass traction. The variables 
have been made dimensionless with respect to a constant reterence 
state Pq, Pq, Tq and Yq. Velocity is referred to the acoustic 
speed Cq, defined by 

c 0 = C/p 0 /P 0 3 1/2 , 

time to t 0 , the homogeneous induction time at the reference state, 
and length to c^tQ. The diffusion terms have been left out because 
they are much too small to play a role in the problem under study. 
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The dimensionless parameters appearing above are the specific— heats 
ratio y, the chemical heat release fi and the activation energy 9. 

Let the reacting gas be confined to the interval O < x < a. At 
the walls the appropriate boundary conditions are 

(2.2) u(0,t) = u (a, t ) = 0. 

The initial state of the gas is taken to be an O(0 -1 ) perturbation 
of the spatially homogeneous and stationary reference state, i.e. , 

(2.3a) u (x ,0) * 0 _1 u 1 (x,O), 

(2.3b) $(x,0) = 1 + e _1 # 1 (x,0) for $ = T, p, Y and p, 

where the precise specification of u^(x,0) and f^(x,0) must await 
the next section. Note that 

(2.3c) Pj(x f O> — p}(x,0) — Tj(x,0) 

in accordance with the gas law (2. le) . An asymptotic solution of the 
initial -boundary— value problem (2.1)— (2.3) is sought in the limit 
9 m, with fi and Y fixed and 0(1), until the localized explosion 
has reached completion. The various stages of evolution are 
detailed in the following sections. 

3. The Induction Stage 

The initial conditions (2.3) suggest that, at least initially, 
the state of the gas remains an 0(fl”^) perturbation of the 
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reference state. During this period, referred to as the induction 
stage, one therefore seeks the expansions 

(3.1) u ~ 9 1 u 1 + ..., f -v 1 + + ..., 

for & = T, p, Y and P 

which, upon substitution into the set (2.1) yield the leading-order 
disturbance equations 

(3.2a) (3/3t ± 3/3x)(p i ± 7 Ul ) * /expCr*), 

(3.2b) a/atcyTjL - (y-i> Pl ] - /exp(T 1 ), 

(3.2c) p t - pi - T j , 3Y 1 /3t = -(l/fl)exp(T 1 ) . 

Except for the nonlinear source term, eqns. (3.2a,b) are simply 
those of linearized acoustics in a uniform atmosphere. It is a 
simple matter to integrate them along the characteristics, as was 
done in [43 and [5] for a different configuration. During 
induction it is enough to concentrate on the variables Tj, pj and 
Uj , because once they are known, the first eqn. in (3.2c) yields Pj 
while the second, combined with (3.2b) and integrated, determines 
Yj according to the expression 

(3. 2d) - [(/-D/npj + flYj = T 1 (x,0) - C ( Y-\ ) /Vlp 1 (x ,0) 

+ <3Y 1 (x,0). 

Equations (3.2) need to be solved numerically, and this was done 


for 
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y = 1.4 


and -for a variety o-f smooth initial conditions and interval lengths 
a. A high-resolution, adaptive ODE integrator was employed to 
integrate along the characteristics. All computations displayed 
thermal runaway, characterized by the unboundedness of Tj and pj 
somewhere in the interval CO,al at a finite time t e . The numerical 
results can all be summarized by considering two representati ve 
cases, for which the initial values of pressure and mass fraction 
correspond to those at the reference state and the initial velocity 
is zero, i . e. , 

(3.3a) p 1 <x,0) = Y & (m v O> = uj(x,0> * 0, 

while the initial temperature perturbations are prescribed as 

(3.3b) Tj(x,0) * bCl-(x/a)l for case I, bCl-(x/a)^3 for case II. 

(Numerical results to be presented below correspond to a *• 0.8, b = 
0.5.) In both cases the initial disturbance has a single maximum 
at x * 0, causing it to become the site of thermal runaway. The 
essential difference between the two cases is that in I the 
temperature disturbance has a nonzero spatial gradient (sharp peak) 
and in II a zero spatial gradient (rounded peak), at x = O. Thus I 
typifies a hot spot located at the boundary (e.g., the shock 
configuration discussed in [43 and C53), and II an internal hot 
spot (easily visualized by a symmetric reflection about the 
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origin). Henceforth the two cases will be referred to, 
respectively, as Type B (boundary) and Type I (internal). Their 
spatial structures, it turns out, are different. 

In the following sections the Type-B problem is discussed in 
detail. The Type-I problem can be treated analogously and is, in 
fact, slightly simpler to analyze; it was deemed sufficient, 
therefore, to simply state its solution in section 6. 

We start with Figure 2, which displays the numerical results for 
the Type-B induction solution. The four graphs there exhibit, 
respectively, the profiles of Tj, p^, Uj and Pj against x for 
increasing values of t, upto the time beyond which the integration 
routine was unsuccessful for a time step 10~^, thus signalling the 
imminence of blowup. An examination of the Tj- profile near blowup 
reveals the birth of a boundary layer at x - 0. Additional 
information is provided by Figure 3, where the function expC- 
Tj(0,t)3 is graphed near blowup. The straight-line graph in the 
figure has slope 1.4 («y), and a t-intercept equal to the blowup 
time t B , allowing one to conclude that 

(3.4) Tj(0,t) -v -jfnCy(t e -t) 1 + o(l) as t •*> t B . 

Figure 4 displays time plots of the solution at x = O, and shows 
clearly that while Tj(0,t) and pj (0,t) become unbounded, Pj(0,t) 
does not. Therefore, p 1 (0,t) must have precisely the same leading- 
order behavior as Tj(0,t), i.e., 

(3.5) p x (0,t) -v -fnCt e -tl + 0(1) as t -*• t e . 


- 10 - 


To summarize, the induction stage exhibits the classic logarithmic 
singularity of spatially homogeneous thermal runaway C83. 

4. Type-B Blowup Structure 

Although numerics has elucidated the temporal character of the 
blowup singularity, further analysis is needed to ascertain its 
spatial structure. This will be done by examining separately the 
boundary layer, whose emergence has already been noted, and the 
region outside. First, it is convenient to introduce a new time 
variable r via the expression 

(4.1) x = t e -t, r > 0. 

Then, following simple manipulations, eons. (3.2a,b) transform into 
(4.2a) dp^/dx - Yd Uj/3r = 0, 

(4. 2b > Y9u 1 /dx - 3p t /3r = (y-l)3p 1 /3r - Yd Tj/Sr *= Xexp(T 1 ), 

where the dependent variables are now treated as functions of x and 

r. The relevant boundary condition is the first of (2.2), re- 
written as 

(4.3) u 1 (0,r) = 0. 

Elementary manipulations on (4.2a,b) and (4.3) yield the following 
integral, which will prove to be of value later on: 
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(4.4) 3>(r) = exp[J exp (Tj (O, t ) ) dt] , 

T 

where 3<r> denotes the disturbance temperature gradient at x * 0, 

Anri 17 ^ 1 tn initial value, i.e.. 

V 

(4.5) 3(0 * [aTj (x,r)/3x] x=0 , 3 Q - 3(t.-) . 

Recall, -from (3.3b) , that vanishes for type I but is negative 
for type B. Then (4.4) shows that J?(r) s O (rounded peak) for the 
former and decreases monotonocal ly to -e (sharp peak approaching a 
cusp) for the latter as r -*■ 0+. 

4.1 The Boundary Layer 

Turning now to the asymptotic analysis near blowup, eqns. (4.2a,b) 
govern the region outside the boundary layer, where the outer limit 
process 


x > O and fixed, r -*■ 0 

applies. The boundary layer, on the other hand, corresponds to the 
i nner limit process 

s > O and fixed, r 0, 

where s(x,r) is the spatial coordinate in the boundary layer, 
reflecting its self— similar structure. The shrinking nature of the 
layer requires x to vanish under the inner limiting process, and 



- 12 - 


then a moment's reflection suggests the definition 
(4.6) s » x/r, 

which assigns coequal importance to the x- and r-derivatives, 
thereby providing the richest equations for the inner limit. (It 
will transpire that this scaling does not quite cover the entire 
boundary layer, but more about that later.) These equations, 
obtained by transforming (4.2a,b) to the (s,r) variables, are 

<4. 7a) (sd/3s - r3/3r) DTj - (/-l)p J 3 = /rexptTj), 

(4.7b) (s3/3s - r3/3r)p 1 + /duj/ds = yrexp(T 1 ), 

(4.7c) (s3/3s - r3/3r)uj + (l/y)3p 1 /3s *» 0. 

It is convenient to isolate the temporal singularity from the 
spatial structure, by setting 

(4.8a) Tj = -jfcn(yr) + f(s,r), 

(4.8b) Pj = -fn(B^r) + g(s,r), 

(4.8c) Uj = h(s,r), 

where the yet unknown constant represents a weak influence of 
the initial conditions on the self-similar boundary layer, and will 
be determined in due course by matching. The structure functions 
f, g and h are assumed to be o(l) in the limit r ■* 0. Substitution 
of (4.8) into (4.7) yields the structure equations 
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(4.9a) 

sf s " 

z * v + 

(/-i)h s 

= e 

(4.9b) 

s(f s - 

g B > - 

r(f r - 

9r> 

(4.9c) 

sh s - 

rh r + 

(1/X)g s 

= 0 


The only boundary condition appropriate -for the above set is the 
Mall condition 

(4. 10) h(0,r) * 0. 

In addition, since the initial data are smooth, the structure 
■functions and their s-deri vati ves are required to be regular in s. 
Consider the asymptotic expansions 

(4.11) § ~ cr 1 (r)$ 1 (s) + or 2 (r> ^2 (s) + ••* * or § “ f , 9, end h, 

as r 0. The gauge sequence <o n (r)> is not yet specified, but a 
clue as to its identity is provided by the integral relation (4.4), 
rewritten as 

X 

(4.12) f s (0, c) * 5 0 r exp[ -J *1/ ( /tlexpff (0, t ) >dt] 

in view of the scaling (4.6) and the prescription (4.8). For small 
r the f-expansion in (4.11) allows the above relation to be reduced 
further to the asymptotic form 

(4. 13) CT t f t ' (0) + ... = 

T 

3 o c X exp[0(l) - {f 1 (0)/X>J* t _1 a 1 (t) dt +...], 


where 
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(4. 14) x = <y-i>/y. 

Recall, from (4.5), that the constant 3q is nonzero for the Type-B 
problem. Then, the assumption that f j ' (0) is nonvanishing 
(involving no loss of generality) leads to the conclusion 

(4.15) or 1 (r) = r X 

if the two sides of (4. 13) are to balance at leading order. With 
determined, it can be shown that the expansions (4.11) proceed 
in powers of r x . 

The boundary-layer analysis can now be carried out, and as 
hinted earlier, the layer is found to have a two-sublayer 
structure. It is convenient to refer to Figure 5 in which the 
various spatial regimes near and beyond blowup are displayed 
schematically. OR refers to the outer region and BL to the 
boundary layer; the latter is subdivided further into an interior 
sublayer Lj and an exterior sublayer (_£. We shall first examine 
the interior sublayer, show that it becomes nonuniform for large s, 
determine the appropriate scaling and expansions for the exterior 
sublayer, and demonstrate that the latter merges smoothly into the 
outer region. Only one or two terms of the expansions in each 
region will be computed; continuation to higher orders is 
straightf orward though increasingly complex algebraically. 

4.1.1 The Interior Sublayer Lj 

Substitution of (4.11) into (4.9) yields the leading— order 
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structure equations for the inner sublayer, 

(4.16a) sfj' - (X+Dfj + ty-Dhj' « 0, 

(4.16b) s(f 1 '-g 1 ') - X(f 1 -g 1 ) - h = 0, 

(4.16c) shj ' - Ahj + (l/y>g 1 * * O. 

The boundary condition 

(4. 16d) hj (0) = 0 

comes from (4.10), and the solution is restricted additionally by the 
requirement that it be regular. If g^ and hj are eliminated from 
(4. 16a— c) , the result is the third-order equation 

(4.17) C(s 3 -s)f 1 'l" + AC (l-3s 2 )f 1 'I * - (s 2 -l/y)f 1 '' 

+ (A+l) (3A-4)sf ! ' + (X 2 -l) (2-A)f ! ■ 0 

for fj. The points s=0 and s=l are singular points of this 
equation and the three linearly independent solutions have the 
asymptotic behavior 

1, s and s£ns as s +0, and 

1, 1 — s and |l-sl 3<y-1)/2y as s •> 1. 

In general one can expect a one-parameter family of regular 
solutions to exist, and numerical computations verify that such is 
indeed the case. A convenient parameter is 
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(4.18) A! = f t (0). 

With -f j known, hj' can be eliminated -from (4. 16a, b) to obtain a 
•first-order dif -f erential equation for gj whose regular solution 
turns out to be 

T 

(4.19a) - C ( X+l/X) / ( X-1 > If j - Cs X /X( Y-l ) 3 J* x -X f 1 *(x)dx, 

0 

and then, (4.16c) integrates to give 

T 

(4.19b) hj = g 1 '/(/X) - Cs X / (XX) ij* x~ X g l "(x)dx, 

0 

where regularity has been imposed again. Thus the full solution at 
this order depends on the single parameter Aj. Graphs of f j, gj 
and hj for Aj « 1 are drawn in Figure 6. 

At this stage the solution (4.B) has the following expansions 
in the interior sublayer: 

(4.20a) *n(Xr) + r X f t (s) + ...» 

( 4 . 20b ) pj ~ -Jn(Bjr) + r x g t (s) + , 

(4.20c) Uj ~ r^j^s) + ... . 

In order to determine the spatial extent of Lj one needs the 
asymptotic behavior of fj, gj and hj for large s. This is easily 
obtained from (4. 17) and (4. 19) , as 

(4.21a) f t ~ -A 1 «Cs (2y-1)/y + s _1/y (A f lns + B f > 
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+ C f s- (y+ l>/y + ... 3 , 

(4.21b) gj a, -A 1 «Cs <2y-1)/y + D g s <y_1)/y + s -1/y (A g jlns + B g > 

+ C g9 - <y+ l>^ + ...3, 

(4.21c> h. - -A 1 «Cs (> '" 1)/y (A hJ en 9 + B«_ ) + C h B _1/y + ...3. 


Here a, Br and C* are constants with values 


a = 2.660, B f * -0.236, C f = 0.0737 

obtained by integrating the equation (4.17) numerically. The 
remaining constants appearing above are given by 

A f = -(2y-i) (y-i) 2 /(2y 3 > , A g - (y+n A f /(y-n , 

A h - 2A f y/(y-i) 2 , B g * A f /(y-i) + (y+i)B f /(y-n , 

( 4 . 22 ) B h - 2B f y/(y-i) 2 - y(3y-i)A f /(y-i) 3 , 

c g - (2y+i>c^/(2y-2) , c h - -3yc f /(y-i>, 

D g * -3C f ^/(y-l) 2 . 

The range of validity of the expansions (4.20) can now be determined. 
For example, substitution of (4.21a) into (4.20a) suggests that the 
latter becomes nonuniform when 

r A s (2y-l)/y = 0(1)f i „e. , s * 0(r“ w ), 


where 


(4.23) 


U 


(y-i)/( 2 y-i) 


—18— 


and the definition (4. 14) of X has been invoked. Correspondingly, 


x = Q (r y/<2y-i), 


0 ( 1 ) . 


The smallness of x indicates that although one has reached the edge 
of Lj, the outer region is still too far. The need for an exterior 
sublayer is therefore apparent. 


4.1.2 The Exterior Sublayer Lg 

In this sublayer the appropriate variables are g and c, with f 
defined by 


(4.24) 


i = r« s = x/r^ 2 *-^. 


The expressions (4.8) for Tj, pj and uj hold again, provided f, g 
and h are now treated as functions of $ and r. The structure 
equations, obtained from (4.9) by transforming from s to ?, are 


(4.25a) (l-Ai)ffg - rf r + (/-l) h| - e*-l , 

(4.25b) <l-/i) g(f f -g f ) - r(f r -g r ) - z u h § ■ 0, 

(4.25c) (l-/i) fh f - rh r + (l//)r w g f = 0. 


hatching requirements imposed by Lj , obtained by substituting 
(4.21) into (4.20) and then employing (4.24), are 


(4.26a) f -v -Aj<x S* 2 ^ l)/y + g( r 2/i j nr ) f 
g ~ -fl,« t ' 2, - u/r * 0( r «), 


(4. 26b) 
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(4.26c) h -v -Ajoc f (y 1)/y CA h (-u £nt + «nf) + B h 3 

+ 0(c 2tf ) , as f + 0. 


DUi UCU oy LMfSae unw Lp SuiutXutl A a auu^n k Im vhw ^ > 


(4.27a) f ~ F Q (f ) + - - - * 

(4.27b) g ~ e 0 (?) + . . . , 

(4.27c) h *v x u llnx H Q ( f ) + Hj (f ) 3 + 


Substitution into (4.25) leads to the di-f -ferential equations 


(l-«) fF 0 ' * exp(F 0 )-l, 

F 0 ' - B 0 * » 0, 

( 1— *i) fH Q * - uH 0 - 0, 

(l-#i)SH 1 ' - iM x = H 0 - (1/7) G q * , 


whose solutions, subject to the matching requirements (4.26), are 


(4.28a) F 0 (f) » Gq( $) * - Int 1+A 1 «f <2y_1 } /y 3 , 

(4.28b) H 0 (f) = -A 1 aC(7-l)/y 2 >f <y “ 1)/y , 

(4.28c) Hj(S) = A t a $ < y ~* > { (27-1 ) /7 2 3 tjfn# + F Q (f)> “ B h 3 * 


r 


i 


Thus the L^-solution can be written as 


(4.29a) Tj ~ — jin ( 7r) + F Q (g> + •**» 

(4.29b) p* -v —In (Bj c) + F Q (f) + ...» 

(4.29c) u t -v tint H Q (f) + H t (f)3 + ... . 
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One must consider the behavior of this solution -for large f in 
order to assess the spatial extent of the sublayer Lg. This can be 
done, -for example, by substituting the large-g behavior of (4.28a) 
into (4.29a). The result is the expansion 

T j ~ -C(2y-l)//>j8n(r y/<2> '“ 1) §) - In^a/) + ..., 
as f + •, 

which clearly becomes disordered when f = 0(r ). 

Correspondingly , x = 0(1), indicating that the edge o-f the boundary 
layer has now been reached. The next step is to see i-f the 
boundary layer merges smoothly with the region OR. 

4.2 The Outer Region OR 

In the outer region, where x and r are the proper variables, the 
solution can be expanded as 

(4.30) fj ~ § 10 (x,t e ) + r f 11 (x) + ..., for fj = Tj, pj and uj, 

where the leading terms are the numerically obtained limiting 
values at blowup and the higher — order terms can be computed from 
(4.2a,b) under the outer limit process. It is a straightf or ward 
matter to establish that a match of (4.30) with the Lg-solution 
(4.29a-c) requires the following asymptotic behavior of the outer 
solution at blowup: 

(4.31a) T t - -l (2y-l)//>jenx - Jn^aX) + ..., 
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(4.31b) Pl -v -C(2y-l)/y>jenx - in<A 1 oB 1 ) + ..., 

(4.31c) u x *>• -A t a x <y_1/y) [{(2y-l)/y 2 >Ct(y-l)/y>JnK 

+ £n (A j«) 3+B^J + . . ., as x -*■ O. 

A careful examination of the numerical solution does, indeed, 
confirm this behavior. The constants Aj and B^, the only ones yet 
undetermined, can then be found by comparing the above expansions 
with the numerical solution. The comparison is made at the "edge" 
of the boundary layer, i.e., for (x,r) satisfying r << 1, r yy * 2y— D 
<< x << 1. It should be emphasized that the structure of the 
blowup singularity is influenced by the initial conditions only via 
these constants; otherwise, the solution has a universal, self- 
similar structure. 

4.3 Summery 

The near — blowup analysis is now complete, and can be summarized. 
In the interior sublayer Lj the expansions are 

(4.32a) T -v 1 + B _1 C-#n(yr) + r^fs) + ...3 + , 

(4.32b) p 1 + fl“*t-ln(Bjr) + r^g^(s) + ...1 + , 

(4.32c) u ~ 0 - *CrNij(s) + ...1 + ..., 

where f g^ and hj are defined by (4.17) and (4.19). In the 
exterior sublayer, the solution is 

(4.33a) T ** 1 + B _1 C-£n(yr) + F Q ( f ) + ...3 + ..., 

(4.33b) p ~ 1 + 0“ 1 C- J en(B 1 r) + F Q (#) + ...3 + ..., 
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(4.33c) u 9 -1 Cr w Inr H Q ($) + Hjtf) + ...3 + ..., 
where Fq, Hq and Hj are given by (4.28). In the outer region the 
expansions take the form 

(4.34a) T ~ 1 + 8 -1 CT 10 (x,t e ) + 0(r)3 + ..., 

(4.34b) p ~ i + © - ^Cpjq(x , t e ) + 0(r)3 + ..., 

(4.34c) u 8 “^Cujq(x *t e ) + 0(r)3 + . .., 

where Tjq, Pjq and Ujg are the terminal values of the induction 

solution, determined numerically. 

The remaining variables p and Y can be computed, upto 0(0“*), by 
appealing to the first equation of (3.2c) and (3. 2d). The results 
are 

(4.35a) 

(4.35b) 

in L j , 

(4.36a) 

(4.36b) 

in l_£, and 

(4. 37a> p ** 1 + 0 -1 Cp lo (x,t e >-T lo (x,t e ) + Q(r)3 + ..., 

(4.37b) Y -v 1 + (0/9) -1 C£(y-l)/y>p lo (x,t e )-T lo (x,t e ) 


fi -v 1 + 8 _1 CJn(y/B 1 ) + ...3 + ..., 

Y *v 1 + (86y)“ 1 C«n(b 1 r) - Fq ( f ) + ...3 + ... 


P ~ 1 + 8 _1 Cfn(y/B 1 ) + r X (gj-fj) + ...3 + ...» 

Y ~ 1 + (8/J/) _1 C£n(b 1 r) + r X £ ( /-l ) g t ~yf + ...3 + ... 


+ b ( 1-x/a) + 0(r) 3 + . . 
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in OR. The constant bj appearing in (4.35b) and <4. 36b) is given 
by 

(4.38) bj « expC/b + YtnY - ( X-l > jBnBi 3 . 

Observe that the BL-solutions (4.32), (4.33), (4.35) and (4.36) 

break down when -inr = 0(8), signalling the end of the induction 
stage, and the onset of explosion. In contrast the OR— solutions, 
(4.34) and (4.37), suffer no disordering and in fact, become 
increasingly accurate as c -*• 0. 
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5. The Type-B Explosion Stags 

The nonuni f ormi ty just encountered decrees that -further 
evolution in the boundary layer occur on the new time scale a, 
defined by 

(5.1> r * e~ 0a . 

For a = 0(1) the limit & •* ® corresponds to a time interval of 
exponential brevity; its role in the evolution of thermal 
explosions was first recognized and exploited by Kassoy C83. The 
two sublayers comprising the boundary layer must again be examined 
in turn. In fact, we shall find that as the boundary layers 
continue to shrink, an expanding void, or an intermediate region 
(denoted by IR in Figure 5) , is created between the sublayer and 
the outer region OR; this region begs a separate treatment. 

5.1 The Interior Sublayer Lj 

The spatial coordinate in this region remains s, now written as 
s * x/r s e& a x, 

thereby expressing explicitly the continuous shrinkage of the 
region. In the (s,cr) variables eqns. (2.1) transform into 

(5.2a) 0" 1 P a + sp s + (pu) s ■ 0, 

PC© -1 u a + su g ] + (l//)p B + puu s ■ 0, 


(5.2b) 
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(5.2c) PC0" 1 T 0 + sT # 3 - €(y-l)/y>C« 1 p a ♦ sp s 3 

+ utpT s - {(/-l)/y>p s 3 - W, 

(5. 2d) fit O’ 1 Y a + sY s 3 + upY s = -(1/0)W, 

(5.2e) p * PT, 

where 

(5.3) W * &~ X pYexp£0( 1 — <j— 1/T) 3. 


The boundary condition (2.2) is rewritten as 


(5.4) u(0,cr) = 0. 


At -fixed s the solution must match with the induction zone as cr 
0. To obtain the necessary conditions one applies the "explosion 
limit" a fixed, Q ® to the Lj-solution (4.32), (4.35) and gets 


(5.5a) 
(5.5b) 
(5.5c) 
(5. 5d) 
(5. 5e) 


T ~ 1+c - P -1 In/ + ... + «Cf 1 (s) + ...3, 
p 1+or & ^ j£nBj + • • . + SCQj (s) + . . . 3 , 
u ^ (s) ^ • 1 f 

p ~ 1 + © _1 Jn(y/ B^ + ... + «Cg 1 (s)“f t (s) + ...3, 
Y -v 1 -a/ (fly) + (B<3y) -1 lnb 1 + ... 

+ fi(fly)“ 1 C(y-i)g 1 (s)-yf t (s) + . . . 3 , a^s cy ^ ^) . 


Here, 

6 - S~ 1 b~ 0&X , 


(5.6) 
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where X was defined in (4.14) and bj in (4.38). These conditions 
reveal that spatial variations in the explosion stage appear only 
at the (exponentially small) Q(fi) level, thereby suggesting that 
the solution is spatially uniform to all algebraic orders in 0. In 

other words, the structure of the interior sublayer consists of an 

extremely weak chemi co-acoustic field superimposed over a uniformly 
exploding atmosphere. Accordingly one seeks expansions of the form 

(5.7a) u -v S Uj(s,a) + ..., 

(5.7b) § -v $q (a; 0) + 6 § 1 (s,tr) + ..., for § * T, p, p, and Y, 

with the understanding that the contain all terms of algebraic 
orders. Substitution into (5.2) finds the $q satisfying the 
standard equations of constant-volume thermal explosion CB1, i.e. , 

(5.8a) 3p 0 /3<r ■ 0, p 0 = PqTq, 

(5.8b) (l/y)P 0 3T 0 /9a ■ -/3P 0 3Y Q /3 <r * W Q a P o Y Q expC0 ( 1— a— 1 /T Q ) 3 . 

The solution, subject to the matching conditions (5.5), is 
(5.9a) § 0 S Q0 + 0 -1 § 01 + ..., for $ = T, p, Y and p, 

where 

(5.9b) T 00 - p 0Q =» (1-<7> _1 , Y 00 » (l + y<5-T 00 )/(yfl) , P 00 * 1, 

(5.9c) T 01 - -(l-<r)“ 2 Jn[y(l-<r) 2 Y 00 ], P Q1 - tniY/ B t ), 

P01 * T 01 + P 0 l /(1 - a) » ***01 * CIiKbi/n-T^l- 


(5. 9d) 
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The structure -functions fj satisfy, to leading order, the equations 


(5. 10a) 

(s3/3s-X)/>j + Pq)U|/3s * O, 

(5. 10b) 

/> 0 <sa/3«-x)u. + a/y)3p«/3s = o. 

(5. 10c) 

(sa/as-X)C/» 0 T 1 -f(y-i)/y>p 1 3 = TjWq/Tq 2 , 

(5. lOd) 

P 0 (s3/3s-X)Y 1 * - T 1 W 0 /(flT 0 2 ), 

(5. lOe) 

Pi ~ p o T i ~ T o p i = °» 


where Wq was defined in (5.8b). Replacement of Tq and Wq by their 
leading— order values from (5.9), followed by the use of the trans- 
f ormations 

(5.11) uj * (1 -a) 1/2 u lt s * s(l-a) _l/2 , 
reduces the set (5. 10) to 


(5. 12a) 

(s3/3s-X)T 1 

+ (7-1) 3^/3* = Tj , 

(5. 12b) 

(s3/3s— X) (T 

j-pj) — 3uj/3s = 0, 

(5. 12c) 

AAA. 

(s3/3s-X)Uj 

+ (l/y)3p 1 /3s = o. 

(5. 12d) 

P 1 * ( Pl~ T l 

)/T C0 , 

(5. 12e) 

(s3/3s-X) 

+ d/(/sy)>T 1 = o. 


Eqns. (5. 12a— c) are identical to (4. 16a— c) if Tp pj, Uj and s in 
the former are identified, respectively, with fj, g^, hj and s in 
the latter. Fallowing the arguments of section 4.1.1, therefore, 
one is led to the solution 
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(5.13a) T l = CAj (a) /Aj3f x (s> , 

(5.13b) pj, = CAj (c)/A 1 Dg 1 (s) 

(5.13c) uj * ( 1— a) 1/2 = (1 -ct) 1/2 CA 1 (a) /A 1 3h 1 (s) , 

where Aj is the constant introduced earlier in (4.18). The 

>S 

amplitude function A^(a) is unknown at this stage, and will be 
determined by matching with the exterior sublayer. So far we only 
know its initial value as a result of matching with the induction 
solution (4.32), i.e. , 

(5. 14) Aj(0) = A 4 . 


It is now a simple matter to solve (5. 12d) for Pj , and compute Yj by 
integrating (5. 12e) subject to the regularity requirement. The 
resulting expressions are 


(5.15a) Pj = (1-®) CAj (cr)/A 1 DCg 1 (s)-f j (s) 3, 

(5.15b) Y j = CA t (<r)/(/3A 1 ) 3C<:(y-l)/y>g 1 (s)-f j (s) 3. 

Both the spatially uniform and the spatially-varying components 
of the expansions (5.7) are thus determined at leading orders, 
although the latter involve A^(a) which is still to be found. It 
is worth noting that the spatial structure of the solution is 
essentially the same as it was at induction-stage blowup p the 
scalings (5.11) simply reflect the temporal evolution of the 
acoustic speed. 

As in section 4.1.1, the Lj-solution breaks down for large s. 
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the nonun iformity now occuring (see the expansion (4.21)) at. 

— 0(1). One is then led to the exterior sublayer. 

5.2 Ths Exterior 5ub!-y*r L £ 

Here the proper variables are f and a where f is now related to x 

A 

and s via the expressions 

(5.16) $ « x e 0oY/ <2Y-1) m (i_or)-l/2 Q e ~ a9u , 

and u was defined in (4.23). In the new variables the full equations 
(2.1) read 


(5. 17a) 

0 1 P a + 

(l-w)f P f + e _aCTW (pu) s = O, 

(5. 17b) 

pe&ou iq- 

* u a + (l-ii)fugi + d/y)pg + puu| » 

(5. 17c) 

Pte" 1 t ct 

+ (i-p)§T|i - c(/-i) /y>t© -1 p a + (i 


+ e~ 9au uCpTg - €(y-l)/y>p f 3 = w. 

(5. 17d) 

PC© -1 Y a 

+ (1-*|)$Y ? 3 + e 9au UPY ? = -(l/fl)W, 

(5. 17e) 

* 

II 

a 



O, 

~/J> Spgl 


where W retains the definition (5.3). The solution is subject to 
the following matching conditions imposed by LjS 


(5. 18a) 

T 

~ T 00 + 9 1 

cT 0 i - P#<2y-i)/y 3 

+ “ • • » 

(5. 18b) 

P 

~ T oo + 0-1 

ct 0 i + p 0 i /( i-^ “ 

p§ ( 2 y-i)/y ] + 

(5. 18c) 

u 

~ - e ~ 9atI P$ 

(y " 1)/y K^ + 



0 -1 CA h jenC(l-or) 1/2 f> + B h l] + ..., 
P ~ 1 + 0 -1 P 01 + . . . , 


(5. lBd) 



-30- 


(5. IBe) 


Y 


Y 00 + (00/) 1 [ln(b 1 /y) - t 01 + 

p§ (2y-i>/yj + as * * o. 


Here, 


(5.19a) P(ff> = aAjtcr) (l-a) <2y 1>/<2 y) , 
with 

(5.19b) P(0) = «A i , 

from (5.14). In obtaining the conditions (5.18) we have employed 
the expansions (4.21) and the solution (5.13); the variables with 
double subscripts are the spatially homogeneous -Functions appearing 
in (5.9). It turns out that compliance with these conditions also 
ensures temporal matching with the induction stage. The L^— 
solution is now sought in the -form 

(5.20a) T -v T oq + 0 -1 T 1 (f,oD + ..., 

(5.20b) p *v T 00 + pj(§,<r) + . .., 

(5.20c) u e~ 0ou Cu 0 ($ T o) + 0' 1 u 1 («,a> + ...3, 

( 5 . 20d ) Y *v Y 00 + 0~ l Yi(e,a) + ..., 

( 5 . 20e ) P -w 1 + 0 -1 p Q1 + ... . 

Substitution into (5.17) shows that (5.17a) is satisfied 
identically to 0(© -1 ). At 0(1), (5.17b) reduces to 
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C 1 — * 1 ) f 3u 0 /3f - *iuq = 0 

whose solution subject to the matching requirement <5. 18c) is 

(5.21) Uq - C(y-l)/y 2 >orPf <y-1)// . 

At O(0 -1 ), (5. 17e) yields 

(5.22) Pi - Tj + p 01 /(l-a), 


while (5.17c) reduces to 

( l-<y> — 2 + (l-*i)$3T 1 /3f - {(y-l)//> [(1-a) -2 
+ Cl-^fdpj/Sf] - Y 00 expC ( 1— cr) 2 T ± □ . 

and, in view of (5.22), simplifies further to 


(5.23) (1-a) -2 + (l-/i) 1 3T t /dS = 7Y 00 expC ( 1-a) 3 




Its solution, consistent with the matching condition (5.18a), is 


(5.24) 


T i = T 01 - (1-a) -2 *n[l + ( 1-a) 2 Pt <2y-1 ' /y ] 


With Tj known, (5.22) defines p t . In order to determine Yj 
consider (5. 17d) at O(0 -1 )j it yields 


-(fly ) -1 (1-a) -2 + (l-w)f3Y 1 /3f 


When linearly combined with (5.23) the above equation leads to 


3ETJ + »y Yjl/3* = 0. 


The matching condition (5. 18e) then provides the following expression 
for YjS 


(5.25) 


/3>'Y 1 = In (bj/X) - Tj . 


It now remains to determine uj, and the function P(cx) (or, 

A 

equivalently, A^(ct)). Both are obtainable from (5.17b) which, at 
Q( 0 - *) , reads 

(5.26) (1-Xi) f 3u 1 /3? - tiUj * -3u 0 /3a - (l/y)3p 1 /3f. 

With Uq and pj known (see (5.21), (5.22) and (5.24)), the 
general solution of the above equation can be written as 

<s.27) uj = t( 2 y-i)/> , >s <) ' _1>/> '[K - {<y-i)/y 2 ><oP) 'in* 

where K(ct) is the integration "constant". As f O, Uj has the 
asymptotic behavior 

(5.28) L 4 •*» C(2y-l)/y>f <y_ 1 )/y [c-{(y-i)/y 2 > (aP) ' 
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+ C(2y-l>/y 2 >P31nf + k] + ... , 

which must agree with the Q(0 - *) term in (5.18c). Matching the #n# 
terms yields the differential equation 

CcrP (<t> 1 ■ = P 

whose solution, subject to (5.19b), is 
(5.29a) P = «Aj, 
or, equivalently, 

(5.29b) = A^l-a)-* 2 '" 1 *" 2 ”. 

With P determined, matching of the $— independent terms in (5.28) 
and (5.18c) yields K: 

(5.30) K * (ocAj/y) C (l/ 2 )Jn(l-a) - fy 2 / (2y-l ) >B h l , 

where the constant B^ was defined in eqn. (4.22). The L^- 
solution at the explosion stage is thus complete. 

It is instructive to compare the solutions in the two sublayers. 
In each the background field is that of a spatially homogeneous 
thermal explosion, but the superimposed spatially-varying field is 
quite different, both in amplitude and structure. In Lj the 
spatial component is exponentially small in amplitude but has a 
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chemi co-acoustic character; all disturbances to the background 
homogeneous field are of the same size. In the spatial 
variations in T, Y and p are O(0 - *) , while those in u and P are 
exponentially small, i.e. , the evolution is essentially due to 
constant-volume chemical amplification of a spatially— non- 
uniform field, with gasdynamics playing a very minor role. 

As or increases, T and p increase and Y decreases, in both the 
sublayers. Eventually, p and T peak when Y 00* the leading term 
in Y, vanishes. This happens at (see (5.9) ) 

(5.31a) a - fly/(l+/5/), 

and the peak values are 

(5.31b) T -v l+/9y, p ~ 1 +fiY. 

At the same time, the O(0 -1 ) term (in T, say; see (5.20a), (5.24) 

and (5.9c)) develops a logarithmic singularity, indicating 
breakdown of the solution and the end of the explosion stage. 

5.3 The Outer Region OR 

This region remains essentially stationary, and hence plays no 
role during the explosion stage. For the sake of completeness, we 
give below the asymptotic form of the outer solution as x -*• 0; 
these expressions are determined by combining (4.31) and (4.34): 

t i + e -1 [-c(2y-i)/y>fnx - jen^®/)] + 


(5.32a) 


• • ■ t 
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(5.32b) p ~ 1 + 0 -1 [-C(2y-l)/y>«nx - IntAtO®!)] + ..., 

(5.32c) u e -1 Aj« x (y_1>/y [{(2y-i)/y 2 >ct-(y-i)/y>#nx 

- £n(A 1 «)3 - B h ] + ... . 

Similar expressions can be written for p and Y. The important 
point to note is that this solution is unmatchable with that in L^; 
for example, to leading order, T is 1 in OR and 1/(1— a) in L^. The 
reason is the emergence of the intermediate region IR in Figure 5, 
created by the receding boundary layer. In this region T must 
vary, at leading order, from the outer value 1 to the inner value 
1/(1— <r) . 

5.4 The Intermediate Region IR 

This region, because of its passive character, will only be 
described very briefly. It is governed by the variables <x and X, 
where X is defined by 

(5.33) x = e~ 0X . 

Matching with the neighboring regions is carried out at fixed <r, by 
setting 


X = -Q 1 *nx 


as one approaches the outer region, and 
X * <r(l-p) - 0 ~ * Jnf, 
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as the boundary layer is approached. Therefore the range of X is 
<5.34) 0 < X < ct< 1— p) ; 

recall that u was defined in (4.23). From (5.20) one can easily 
conclude that leading-order matching with L E requires 

<5. 35a) T ** 1/ < l-o) + ... , as X -*■ o<l-p>, 

with analogous expressions for p and Y, while 

(5.35b) p -w 1, u e ~®X<X-l)/y t<y-l) 

Therefore the solution is sought in the form 

(5.36a) * •* f° + ..., for i ■ T, p, P and Y, 


and 


(5.36b) 


.. m -ox(y-i)/y ..0 . 

U 'V G U T . a a « 


In the (X,<r) variables the full equations (2.1) read 
P a - (up) x « o, 

p[u a - uu x e -®<»- x >] - (l/y)p x * 0, 

p[T a - uT x b “® (<t " X> ] - C(y-l)//>[p a - up x e “* <a “ X> ] - W, 

p[Y a - uY x e _0<tT_X) ] - -<l/j9)W t 
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P = PT , 

where W retains the definition (5.3). Clearly, the solution is 
stationary, i.e. , independent of a to all alyteLir aic orders. 
Specifically, the reaction term W is exponentially small since one 
expects T < l/<l-o). The leading-order terms can then be 
determined simply by appealing to the matching conditions (5.35), 
and one finds that 

(5.37a) T° ■ l/Cl-(l-w) -1 X3, 

with analogous expressions for p and Y, while 

(5.37b) P° ■ 1, u° - t(/-l)/y 2 >A 1 a (l-*i) -1 X. 

It is a simple matter to check that the above solution also matches 
with the outer expansions (5.32) as X 0. 

The analysis of the Type-B explosion is thus complete. 

6 . The Type-! Explosion 

In this section details are largely omitted and emphasis is on 
the results, since the treatment follows closely the Type-B 
analysis just concluded. 

6. 1 The Induction 8tage 

Figures 7(a-d) display the numerical solution of the induction 
problem. The graphs are self-explanatory. Similarity with Figures 
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2(a-d) is obvious, but two points of contrast are noteworthy. 

First, the temperature pro-file now has a rounded peak. Second, the 
boundary layer is thicker; this can be seen more clearly in Figure 8, 
where Tj (x ,t) /T j (0,t ) is plotted at the last successful time step 
for each of the two cases. 

The boundary layer retains the form <4. 8) and a two— sublayer 
structure emerges once again. The expansions are 

(6.1a) T -v 1 + 8“ 1 [-ln(yr) - zinz + 

r {-AjS 2 + A 2 > + . ..] + • ••» 

(6.1b) p 1 + 8 -1 [-^(Bjr) - rinr t ( X+l ) /yjAj + 

r C-Ajs 2 + f(y+l)/(y-l)>A 2 + A t /Yl + ...] + ..., 
(6.1c) u ~ 0 _1 [rfnr 1-2 Ajs/yi + 

z C2A 2 /(y-i) - Ai/yis + ...] + ... 

in Lj , and 

(6.2a) T a, 1 + © _1 [-j|n(yr) - Ind+AjC 2 ) + ...] + ..., 

(6.2b) p 1 + 8“ 1 [-ln(B 1 r) - Ind+Aji 2 ) + ...] + ..., 

(6.2c) u ~ 0 _1 [r 1/2 Jnr {-(2/y)AjiC> + r 1/2 C-(2/y)AjC in(l+A t < 2 ) 

+ C2A 2 /(y-l) - Aj/yJCD + ...] + ...» 

in Lg-. The coefficients Aj , A 2 and Bj are to be determined by 
matching with the outer solution as before. The spatial coordinate 
£ in l_ E is defined by 

(6.3) C - x/r 1/2 . 
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implying that the boundary layer is now 0(r 1/2 ) thick, and hence 
thicker than the 0(r ') Type-B layer. The -finding o-f Figure 

8 is thus confirmed. 

It turns out further that the solution is uni fur inly valid all 
the way to ( = 0, so that the interior sublayer is, in fact, 
superfluous. 

For smooth merging with the boundary layer the outer, 
numerically computed solution is required to have the asymptotic 
form 

(6.4a) T -v 1 + 0" 1 C-2*nx - inOT^) + ...3 + ..., 

(6.4b) p *v 1 + 0 ~ * C-2£nx - £n(A^B^) + ...1 + ..., 

(6.4c) u *** 0 -1 C-4(A 1 /y)xlnx + €-2 ( Aj//) jfnA x 

+ 2A 2 /()'-l) - A 1 /y>x + ...3 + ..., as x ■* 0. 

This behavior was confirmed, and the constants Aj, A 2 and Bj 
computed, by comparing the numerical solution with the above 
expansions. The remaining variables and Yj can be determined 
as before, by appealing to the first member of (3.2c), and 
(3. 2d) . 

6.2 The Explosion Stage 

The analysis proceeds as in section 5.2. The appropriate 
coordinates are a and C, and the requirement of matching with Lj 
is replaced by the condition of regularity at ( ■ 0. The solution 


turns out to be 
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( 6 . 5a > 

T 

** T oo + 0 1 | 

[ T oi _<1_a)-2j * ntl+A i* 2> ] + 

( 6 . 5b ) 

P 

- t 00 + e- 1 | 

[ T Oi-(l-er )_2j fn€ 1 +A 1 C 2 >+Poi / <l-^ ) ] + •••» 

(6.5c) 

u 

-v B~ 0a/2 [2rt 1 C<r/fy<l-CT) >] + ... 


(6. 5d) P - 1 + B -1 P 01 + 

(6.5e) Y -v Y 00 + (©fl/)" 1 [ln(b 1 /y)-T 01 

+(l-a) -2 Jn€l+A 1 C 2 >] + ..., 

where the constant bj appearing in (6.5e) was defined in (4.38). 

The doubly subscripted quantities correspond to the spatially 
homogeneous explosion, and were introduced in (5.9). The explosion 
stage peaks just as it did for Type-B, and the remarks at the end 
of section 5.2 remain valid. Finally, the IR— analysis of section 
5.4 carries over, with obvious modifications. 
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7. Concluding Remarks 

The spatial structure and temporal evolution of a localized 
thermal explosion in a confined gas have been described 
mi* thematically. Localization, rather than a spatially uniform 
explosion, occurs as a result of system nonhomogeneities, here 
modelled by a slightly nonuniform initial temperature. Attention 
is confined to what may be called the fast-reaction limit, 
characterized by the initial induction time of the reaction being 
comparable to the initial acoustic time across the vessel, so that 
diffusion plays no role. This limit can be achieved if the initial 
temperature of the unreacted gas has been raised to a sufficiently 
high level, perhaps by the passage of a strong shock. By contrast, 
the slow-reaction limit would correspond to the induction time and 
the conduction time being of the same order. The latter problem 
was the subject of Poland and Kassoy's investigation C63. 

The explosion is shown to develop in two distinct stages. The 
first stage is induction, characterized by small perturbations 
about a spatially uniform state, where the primary interaction is 
between linearized acoustics and weak but nonlinear chemical 
heating. Chemical amplification leads to localized thermal 
runaway, or blowup of the perturbations, at a time and location 
determined by the initial and boundary conditions. The spatial 
structure at blowup is self-similar, differing slightly depending 
upon whether the runaway site is at the boundary or in the interior 
of the domain. 

Induction is followed by explosion, characterized by 0(1) 
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variations in the state o-f the gas. The characteristic chemical 
time plunges dramatically. The acoustic time drops as well, but 
not nearly in the same proportion, so that explosion is dominated 
by chemical heating. There is no time -for expansion, with the 
result that changes in the velocity and density -fields are 
negligible. Thus the gas explodes locally at essentially constant- 
volume conditions, with little change in the spatial structure that 
it inherited at runaway. (Analysis in the Appendix shows that i-f 
thermal expansion is admitted, the corresponding spatial structure 
is necessarily singular.) The explosion stage ends when temperature 
and pressure within the explosion have peaked, the -final values 
being exactly the same, to leading order, as in the spatially 
homogeneous case. The subsequent expansion of the hot, highly 
compressed gas, and the eventual development of a blast wave, are 
currently under study. 
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Appendl* 

The setup (4.8), employed in the text -for analyzing the spatial 
structure of blowup, is based on the numerical observation that 
both Tj and pj exhibit identical, -Inr behavior as r + 0. 

This observation, found to hold for all the numerical runs 
undertaken, implies that blowup is a const ant -volume process, since 
density perturbation = pj - Tj remains bounded. 

Let us now consider the possibility that for some initial 
conditions, blowup lies partway between a constant— volume and a 
constant-pressure process, and ask whether a self -similar structure 
consistent with this notion exists. Accordingly, we replace (4.8) 
by 

(A. la) Tj ~ -inr + f 0 (s) + ..., 

(A. lb) Pj -v -Ainr + gQ(s) + ..., 

(A.lc) Uj ** hQ(s) + ..., 

where 


0 < A < 1. 


The case A » 0 corresponds to a constant-pressure situation, and 
A ■ 1 to the constant-volume case already discussed. Substitution 
into (4.7), followed by some rearrangement , yields the leading- 
order structure equations 
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(A. 2a) s<l-s 2 )f 0 ' + Cl-s 2 -€<y-l)/y>AD - <l-ys 2 )exp f Q , 

(A. 2b) <l-s 2 )g 0 ' - sCA-yexp f 0 3, 

(A. 2c) y<l-s 2 )h 0 ' - yexp f Q - A. 

The transformation 

(A. 3) f Q = -InF 

reduces (A. 2a) to the linear equation 

(A. 4) s ( 1— s 2 ) F ' - Cl-s 2 -€(y-l)/y>A3F * i-ys 2 . 

Once F is known , f 0 , g 0 and hQ can be computed sequentially from 
(A. 3) and <A.2b,c). 

Equation (A. 4) has singular points at s — 0 and 1. It can be 
shown that in general integration can remove at most one 
singularity, thereby yielding solutions which are singular 
either at 0 or at 1 . Such solutions can evolve only from very 
special, singular initial conditions, and are therefore 
unacceptable if the initial data are smooth. The only regular 
solution is the constant 

f - 

which requires 

A = 1, 

corresponding to the constant-volume blowup already discussed. 
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Figure 1 

A scematic of the temperature profiles for the Type-B (sharp 
peaked) and the Type-I (round-peaked) explosion. 
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Figure 2(a) 

Type-B induction-stage profiles for (a) Tj t (b) pp (c) Uj and 
(d) Pj. The profiles are plotted at (i) t=0, (ii) t=0.2, 

(iii) t=0. 4, (iv) t=0.43, (v) t=0.442, (vi> t=0. 44625. Estimate 

of blowup time is t_ = 


0.446890. 
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Figure 2(b) 
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Figure 2(c) 
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Figure 4 


Type— B plots of Tj(0,t)» Pj (0,t) and Pj (0,t) . 
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Figure 5 


A schematic o-f the spatial zones at and beyond blowup. 

BL: Boundary Layer , Lj : Interior Sublayer, LjrS Exterior Sublayer 
OR: Outer Region, IR: Intermediate Region. Not to scale (The t- 
scale is cosiderably stretched) . 
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Figure 6 


Profiles o-f structure -functions -f < . g« and h< for +< (0) 
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Figure 7 (a) 

Type— I induction-stage pro-files for (a) Tj, (b) pj, <c) Uj 
<d) Pj. The profiles are plotted at (i) t=0, (ii) t=0.2, 
(iii ) t=0. 4 , (iv) t=0.422, (v) t=0.432, (vl) t=0.435S75. 

Estimate of blowup time is t e = 


and 


0.435SB0. 
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Figure 7(b) 
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